Restricting profile function of hedgehog Skyrmion 



Jun Yamashita 1 and Minoru Hirayama 2 

1 Department of Physics, University of Toyama, Gofuku 3190, 

Toyama 930-8555, Japan 
Email : fxjun@jodo.sci.u-toyama.ac.jp 

2 Department of Physics, University of Toyama, Gofuku 3190, 

Toyama 930-8555, Japan 
Email : hirayama@sci.u-toyama.ac.jp 

June 2006 

Abstract 

The profile function for the hedgehog Skyrmion is investigated. 
After discussing how the form of the profile function is restricted by 
the field equation, the static energy is numerically calculated. It is 
found that the profile functions considered here sometimes give the 
static energy smaller than previous ones. 

The Skyrme model is an attractive candidate for the effective theory of 
Quantumchromodynamics (QCD)[21 EH- Its static version is defined by the 
static energy functional 

o 



(1) 

R, ={diU)U\ (2) 
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where U = U(x) is an element of SU(2). We are here adopting the nor- 
malization of E used in refs. |U El E] • The important topological quantity 
associated with the Skyrme model is the baryon number 

B = J d 3 x trilURjRk). (3) 

With these definitions, we have the Faddeev-Bogomolny inequality E > \B\. 

We here restrict ourselves to the case of B = 1 and assume that U(x) 
obeys the hedgehog Ansatz: 

U(x) = exp [if(r)x • r] , r = \x\, (4) 

where r denotes the triplet of Pauli matrices. The true profile function f(r) 
minimizes E. It is known that E can be written as [6 

i r°° 

E =— / s(r)dr, (5) 
3vr Jo 

e(r) =r 2 f(r) 2 + 2 sin 2 f(r) [f (r) 2 + 1] + ( 6 ) 



where f'(r) denotes the derivative of f(r) with respect to r. The field equa- 
tion for f{r) is given by 



[r 2 + 2sin 2 /(r)] f"(r) + 2r/'(r) + sin2/(r) 



/» 2 - 1 - ^ 



0. 

(7) 



The static configurations which minimize E for a given B is called Skyrmions. 
In the case of B =1, Battye and SutcliftV4 began their numerical analysis 
with the initial trial configuration 

f(r) = 4Arctan(exp[— r]) (8) 

which gives 22=1.24035. It was concluded that the energy of the B = 1 
Skyrmion is equal to 1.232[4j. In the later literatures, we find 1.2322 
E]. The value 1.2322 is sometimes referred to as "true" [7J E]. Krusch[9 
investigated Skyrmions on the 3-sphere 5* 3 with the radius L. In the flat 
space limit (L —>■ oo limit), he obtained the profile function 

f(r) =n - 2Arctan(/cr) (9) 
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where k is an adjustable constant. 



In this paper, we first discuss how the functional form of the profile func- 
tion is restricted by the field equation. It is then important to know how the 
adjustable parameters come in the profile function. We investigate a method 
to introduce the two parameters in the profile function. We then calculate 
numerically the minimal value of the static energy under the variation of 
the two parameters. Though numerical calculations aided by Mathematica 
package, we obtain the static energy of hedgehog Skyrmion comparable with 
or smaller than the above values. 



To see the mathematical properties of f(r) implied by the field equation 
involving a transcendental function of f(r), it is convenient to transform 
it into an algebraic differential equation. It is also convenient to define an 
independent variable which varies in a finite interval as r ranges from to 
oo. Thus we introduce v(z) and z by 



v(z) =tan f(r), 
z 



Then we have the field equation 



d 2 v _ 1 
dzi~2 
1 
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v — 1 V v — z 

1 1 2 
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z v — z 



r 2 + 2' 



fdvV 
\dz) 

dv v [v(z + 1) - 2z] 
dz 2z 2 (z - if (z - v) 



(10) 
(11) 



(12) 



and the static energy 



E 



(z — v) \J z (1 — z) f dv\ 



v (v — 1} 



dz J ^ 



v (3zv — 4z + v ) 



2(v-l) 



z 3 (1 - z) 



dz. 



(13) 



The field equation Eq.(|12|) allows us to assume that, in the neighborhood 
of Zq, there is a solution of the form 



v(z) = /J Vj(z — z ) a+j , Vj : const. (j = 0, 1, 2, 

j=0 



(14) 
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Then the procedure known as the leading order analysis in the theory of 
nonlinear differential equation is applicable ^Uj- If we assume a < and 
z ^ 0, 1 and suitable the above expression for v(z), the l.h.s. of Eq. (fT2"j) 
becomes —a (a + 2) t> (z — zo) a ~ 2 /2 + • ■ • , where • • ■ represents terms less 
singular at z = zq. We see that a must be —2. Though similar considerations 
in the cases z = and z = 1, we conclude 

a = 1 if zo = 0, 

a = 2 if «o = 1, (15) 
a = -2 if z j£ 0, 1. 

We also find that, by Eq.(fT2"j). v i, t>2, t>3, • • • are determined successively in 
terms of z and v . Note that the behavior of v(z) implied by Eq.(JT3J) matches 
the boundary condition 

/(0)=7T, /(oo) = (16) 

since z = 0, 1 correspond to r = 0, oo, respectively. With the aid of the 
boundary condition (JTfiJl and the continuity of f(r), we see that there should 
be z such that 



/(ro)=f, ^0 = ^^, <«,<!■ (17) 
Noting that Eq. p7jl implies v(zq) = oo, we find that v(z) takes the form 

where w(^) is a smooth function which is nonvanishing and finite for ^ 
z ^ 1. Then we see that the natural relation between /(r) and v(z) inverse 
to Eq. lfTUJ) is given by 

f vr - Arctany/^z), ^ 2 ^ z , , 1fV . 
I ArctanA/f (z), z S z S 1. 

Although it is difficult to obtain an exact u>(z), we know that it must be a 
moderate function for ^ z ^ 1. Thus we may be allowed to truncate W (2:). 
If we set w(z) as 

JV 

to(z) = ^2w k {z - z ) k , JV:small, (20) 

fc=0 
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and still use the field equation for w(z), then wi,w 2 ,--- are successively 
determined as 

(4z - 1) w 

Wi = ; r, 21) 

(148z 2 - 72z + 25) w + 16 (3 - z ) 

W2 ~ AO 2/1 \2 ' 

A&z 1 (1 - zq) 

(51 - 20(k + 292z 2 - 4084) + 32z (3 - 9z + 2z 2 ) 
W3 = 96^-1)' ' <23) 

and so on. In this case, v(z) involves two adjustable parameters zo and wo- 
We report here the results for some small N. We first consider the simplest 
case that w(z) is constant, i.e., N = 0: 

w(z) = w . (24) 

The behavior of the profile function near r = 0, r , 1 is calculated to be 

/(r)~7r-Kir, «i = — W^, (r~0), (25) 

V 2 



f( r ) ~ ?-«2(r-r ), K2 = \ — — — , (r~r ), (26) 
2 y w 

/(r) 2 , k 3 = , (r ~ oo). (27) 

We have numerically estimated E with changing wq and by the step 0.001. 
We find that 

E = 1.23186- •• = 1.2319 (28) 

for 

. . 0.673^(1 -zf 
v ( z ) = i > 29 

V ; (2- 0.279) 2 V ; 

is at least a local minimum. The profile function and energy density of this 
case is depicted in Fig.l and Fig. 2, respectively. 
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Fig.l: Profile function f(r) for 
w = 0.673, z = 0.279. 




Fig. 2: Energy density e(r) for 
w = 0.673, z = 0.279. 

In the N — 1, 2, 3 cases, with the help of (j21H23|) . we have 



N = l, 
N = 2, 
iV = 3, 



w 
w 
w 



0.670, 
0.485, 
0.385, 



0.275, 
0.327, 
0.249, 



E 
E 
E 



1.23215- •• 
1.34000 • • • 
1.34234- •• 



1.2322. 
1.3400. (30) 
1.3423. 



We see that the N = case gives a smaller E than N = 1, 2, 3 cases. 

On the other hand, if we do not make use of (j21H23|) and consider the case 

that w(z) is a polynomial of z of second order, we have the minimal energy 



E = 1.23147- •• ^ 1.2315 



(31) 
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for 



( 1 1^ 

ute) = Z{ ~ Z> (0.625 + 0.3222 - 0.356z 2 ) 
^ - 0.278) 2 V ; 



(32) 



If we add some terms of the form of z k , k ^ 3 to the above (0.625 + 0.322,2 — 0.356,z 2 ) 
we would obtain smaller E's. We have thus suggested that the value of the 
static energy may become smaller than 1.2322. 

In ref.[2], some physical properties of nucleons were discussed in terms of 
the numerical solution of Eq. (J7J) and the input data that masses of nucleons 
and delta particles are given by M N = 939 MeV and M A = 1232 MeV, 
respectively. We briefly compare here the results obtained in the above N = 
case (v (z) of Eq. (}2*9~j) ) with those in ref.[2 . In ref.[2], Mj^ and Ma are given 
as M N = M + 3/(8A), M A = M + 15/(8A), where the parameters M and 
A are related to E of Eq.© and /(r) by 



M 
A = 



3n 2 F w 



-E, 



e 
2tt 

3e 3 F w ' 



A = 8 / dr r sin f(r) 
Jo 



1 + 



df(r) 
dr 



+ 



sin 2 /(r) 



(33) 
(34) 

(35) 



where e and F n are coupling parameters appearing in ref. [2]. The input data 
for M N and M A yield M = 866MeV and A = 0.00514(MeV)" 1 up to three 
decimal places. The isoscalar electric mean square radius (r 2 )j=o and the 
isoscalar magnetic mean square radius a/ (^ 2 )m,/=o are defined by 



dr r 2 sin 2 f(r)f 



(r 2 ) 



M,i=0 



2 dr r* sin 2 f(r)f'{ 



eF n y J °° <ir r 2 sin 2 f(r)f'(r) 



(36) 
(37) 



In ref.[2], the results of the analysis are stated in the following way. 



M = 36.5—, A = 50.9, e = 5.45, = 129 MeV, 
e 



V(r 2 )i=o = 0.59 fm, \ (r 2 ) MJ=0 = 0.92 fm 



(38) 
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while we have in the N = case 

M = 36.5—, A = 52.2, e = 5.48, F„ = 130 MeV, 

e 

vV>/=o = 0.586 fm, yj (r 2 ) MJ=0 = 0.920 fm. (39) 

We find that our AT = case almost reproduces the results of ref. 0. The 
characteristic feature of our procedure is that we make use of v(z) consisting 
of two branches of the function arctan ^v(z). It does not cause any difficulty 
in the behavior of the profile function f(r) as is seen in Eqs. (j25H2*7j) . Because 
of the intimate relationship between the Faddeev and Skyrme models [TTj . 
the minimal energy of the soliton of the Faddeev model ^2] with Hopf charge 
1 may be made smaller. 
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